C*********************************************************************C
C*                                                                   *C
C*  chamang.for                                                      *C
C*                                                                   *C
C*  Written by:  David L. Huestis, Molecular Physics Laboratory      *C
C*                                                                   *C
C*  Copyright (c) 2000  SRI International                            *C
C*  All Rights Reserved                                              *C
C*                                                                   *C
C*  This software is provided on an as is basis; without any         *C
C*  warranty; without the implied warranty of merchantability or     *C
C*  fitness for a particular purpose.                                *C
C*                                                                   *C
C*********************************************************************C
C*
C*	To graph the accuracy of Chapman Function approximations vs
C*	solar zenith angle for X = 800
C*
C*  USAGE:
C*
C*	FORTRAN chapang.for
C*	FORTRAN chapman.for
C*	FORTRAN chaplit.for
C*	LINK chapang+chapman+chaplit
C*	chapang > chapang.out
C*	GRA .lgo chapang.kvs
C*
C*  EDIT HISTORY:
C*
C*	01/29/2000 DLH	Released version
C*
C**********************************************************************

	implicit real*8(a-h,o-z)
	real*4 atm_chap_num
	parameter (nc=10)
	dimension calc(nc), err(nc)
	common/atm8_chap_cm/Fn(0:3)
	X = 800
	  DO j=1,90,1
	    chi0 = j
	    temp = atm8_chap_asy(X,chi0)
	    calc(1) = Fn(0)
	    calc(2) = Fn(3)
	    calc(3) = atm8_chap_deq(X,chi0)
	    calc(4) = atm_chap_num(sngl(X),sngl(chi0))
	    calc(5) = atm8_chap_Ch31b_22( X, chi0 )
	    temp = atm8_chap_Ch31b_38( X, chi0 )
	    calc(6) = Fn(0)
	    calc(7) = Fn(3)
	    calc(8) = atm8_chap_Fi64_3(X,chi0)
	    calc(9) = atm8_chap_GLG64_11(X,chi0)
	    calc(10) = atm8_chap_Sw64_52(X,chi0)
	    truth = atm8_chap_num(X,chi0)
	    do k=1,nc
	      err(k) = abs(truth-calc(k))/truth
	    end do
	    write(*,2400)j,(err(k),k=1,nc)
	  end do
2400	format(i6,1p6d10.3/(6x,1p6d10.3))
	end